Effectiveness and improvement of cylindrical 
cloaking with the SHS lining 

Allan Greenleaf 
Yaroslav Kurylev 

Matti Lassas 
Gunther Uhlmann 

Abstract 

We analyze, both analytically and numerically, the effectiveness 
of cloaking an infinite cylinder from observations by electromagnetic 
waves in three dimensions. We show that, as truncated approxima- 
tions of the ideal permittivity and permeability tensors tend towards 
the singular ideal cloaking fields, so that the anisotropy ratio tends to 
infinity, the D and B fields blow up near the cloaking surface. Since 
the metamaterials used to implement cloaking are based on effective 
medium theory, the resulting large variation in D and B will pose a 
challenge to the suitability of the field averaged characterization of e 
and fi. We also consider cloaking with and without the SHS (soft- 
and-hard surface) lining, shown in [6j to be theoretically necessary 
for cloaking in the cylindrical geometry. We demonstrate numerically 
that cloaking is significantly improved by the SHS lining, with both 
the far field of the scattered wave significantly reduced and the blow 
up of D and B prevented. 

1 Introduction 



1.1 Background and history 

There has recently been much activity concerning cloaking, or rendering ob- 
jects invisible to detection by electromagnetic (EM) waves. For theoretical 
descriptions of EM material parameters of the general type considered here, 
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see [U El El lU El E] ; for numerical and experimental results, see [3 El El QUI EH] • 
Related results concerning elastic waves are in [121 EH EH] • All of these pa- 
pers treat cloaking in the frequency domain, using time harmonic waves of 
some fixed frequency k > 0; this is not unreasonable, since the metamaterials 
used to implement these designs seem to be inherently prone to dispersion, 
for both practical and theoretical reasons [T51 ESI II]- See [TT| for a treat- 
ment of cloaking in the time domain. One can also design electromagnetic 
wormholes, which allow the passage of waves between possibly distant points 
while most of the wormhole remains invisible [TH [T9] . 

When physically constructing a cloaking (or wormhole) device, one is of 
course not able to exactly match the ideal description of the EM material 
parameters (electric permittivity e and magnetic permeability fi, for the pur- 
poses of this paper). Any actual implementation will only realize a discrete 
sampling of the values of e and /i, and not be able to assume the ideal values 
at points x on the cloaking surface(s), where the tensors e(x) or (i(x) have 
or oo as eigenvalues. 

1.2 Approximate cloaking and linings 

The purpose of the current paper is twofold. First, we wish to explore the 
degradation of cloaking that occurs when the ideal material parameter fields 
are replaced with approximations obtained by limiting the anisotropy ratio, 
L, as described below. This was studied in two very interesting recent papers. 
Ruan, Yan, Neff and Qiu [20] consider the effect on cloaking of truncation of 
the the material parameter fields, while Yan, Ruan and Qiu [21], study the 
effect of using the simplified material parameters employed in [HI E] . In [20J , 
it is shown that cloaking of passive objects, i.e., those with internal current 
J = 0, holds in the limit as L — > oo, but a slow rate of convergence of the 
fields is noted. The current paper reproves this and demonstrates the blow 
up of the B and D fields at the cloaking surface as L — > oo. 

Secondly, we consider the effect of either including or not including a 
physical lining to implement the soft-and-hard surface (SHS) boundary con- 
dition, which is a boundary condition originally introduced in antenna design 
[221 1231 l2lj . As we proved in |6J, cloaking EM active objects, i.e., objects 
with generic J ^ 0, imposes certain hidden boundary conditions on the waves 
propagating within the cloaked region. Any locally finite energy wave satis- 
fying Maxwell's equations in the classical, or even weak, sense must satisfy 
these conditions. We note that in our terminology, the fields (E, H, D, B) are 
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a finite energy solution if all the components Ej, Dj, Hj, and Bj are locally 
integrable functions; the energy of the fields is locally finite; and they satisfy 
Maxwell's equations in the classical or weak (distributional) sense. The rea- 
son why we concentrate on such solutions is that the effective medium theory 
of metamaterial requires that the scale at which the EM fields change signif- 
icantly is larger than the size of the components (or cells ) used implement 
the metamaterial. 

For a cylinder cloaked by what is called the single coating construction 
in [6], and which corresponds most closely with the cloaking considered in 
[H [H [8] , the hidden boundary conditions are the vanishing of the angular 
components of E and H . This is exactly the SHS condition associated with 
the angular vector field J|. We show that using a SHS lining has two benefits: 
blow up of B of the cloaking surface, which may seriously compromize effec- 
tive medium theory for metamaterials, is prevented and secondly the farfield 
pattern of the scattered wave is greatly reduced. It is shown in [6] that there 
is no theoretical, frequency-dependent obstruction to cloaking, but with cur- 
rent technology, cloaking should be considered as essentially monochromatic, 
and we will work at fixed frequency k. 

2 Single coating of a cylinder 

Let us consider Maxwell's equations on M 3 , 

V x E = ikB, 

V x H = -ikD, 
D = eE, 

B = fiH, 

where for simplicity we have taken the conductivity a = 0. 

We consider here EM waves propagating in metamaterials, which allow 
one to specify e and \i fairly arbitrarily. These are typically assembled from 
components whose size is somewhat smaller than the wavelength. Ideal mod- 
els of cloaking constructions consist of prescribed ideal parameter tensors 
e, fi, describing coatings making objects invisible to detection by waves of 
frequency k; physically, these would be implemented using metamaterials 
designed to have e,fi as effective parameters (at the specified frequency). 
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Note that in the cloaking constructions the ideal parameters are singular on 
a surface surrounding the object, the cloaking surface. As a result, discussed 
further below, we need to consider Maxwell's equations holding not only in 
the classical sense but also in the sense of Schwartz distributions |25j . 

In the following, we describe the non-existence results for finite energy 
distributional solutions with respect to the ideal parameter fields, the con- 
sequences of approximative material configurations, and the role of the SHS 
lining. 



2.1 Equations for an ideal single coating 

On M 3 , with standard coordinates x = (xi,X2,x^), we use cylindrical co- 
ordinates (r, 9, z), defined by (r,8,z) i— > (r cosd, r sin 6, z) 6 IR 3 . In [B] we 
considered Maxwell's equations on M 3 \ S, 

Vx£ = ikB , V x H = -ikD + J, 
D = eE, B = jJH, 

where e and Jl correspond to the invisibility coating materials on the exterior 
of the infinite cylinder N 2 = {r < 1} and are Euclidian inside N 2 . e and fx 
are singular at S, namely, for r = | (xi, x 2 ) | — > 1 + , 

max M4 = ((r-l)- 2 ) ^oo, 

where \j(x),j = 1,2,3, are the eigenvalues of s(x) or In particular, 

we considered the question of when there are fields E, H, D, B that together 
constitute a finite energy solution of Maxwell's equations in the sense of 
distributions. It was shown in [B] that, in the presence of internal currents J 
when the cloaked region is, e.g., a ball, such solutions do not generally exist. 
Let us discuss why this is so. Even for cloaking passive objects, i.e., J = in 
the cloaked region, the singular material parameters give rise to solutions in 
Maxwell's equations that correspond either to surface currents (see below) 
or to the blow up in the fields at the cloaking surface. Thus, if the material 
does not allow such currents to appear, then the resulting fields must not 
blow up. 

Let us next consider in the scattering of a plane wave by a cloaked cylin- 
der, that is, the case when we have no internal currents and the EM fields 
have asymptotics at infinity corresponding to a sum of a given incident plane 
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wave (E m , H m ) and scattered wave (E sc , H sc ) that satisfies the Silver- Miiller 
radiation condition [2B]. It was shown in [6] that with respect to cylindrical 
coordinates (r,6,z), 

lim eg(x) • E(x) = 0, lim eg(x) • H(x) = 0, 

r—>l+ r— +1+ 

where eg is the angular unit vector. Let e z be the vertical unit vector. For 
general incoming waves, we have that 

lim e z (x)-E(x) - a e (^) = 0, (1) 

r-+l+ \x\ 
~ X 

lim eJx) ■ H(x) — dh(-, — r ) = 

where a e and ah do not vanish. In the treatment of cloaking passive objects 
[H E] it is assumed a priori, based on the behavior of rays on the exterior, 
that the inside of the cloaked region is "dark", that is, the fields E and H 
vanish in {r < 1}. (However, see also [201 EI], where the behavior of the 
fields within the cloaked region is studied.) Under this assumption, the E 
and H fields have jumps across E, 

b e = (v x E^j | s+ - (i> x E^j | s - = a e (x)e e , 

b h = (v x H^j | s + - (v x H^j | s - = a h (x)e e . 

(Here v is the Euclidian normal vector of S, which is just the radial unit 
vector e r .) This implies that 

V x E = ikB + K sur j, V x H = —ikD + J sur f, 

in the sense of distributions on R 3 , where the singular distributions K sur f = 
b e ^E, Jsurf = are terms that can be considered either as magnetic and 
electric currents supported on E, or, as below, idealizing the blow up of D 
and B near E. Here, 5s is the distribution defined by 

f(x)5?.dx= \ f(x)dS(x), 
is 

where dS is the Euclidian surface element on the surface E, for any smooth 
test function /. We refer to such strongly singular field components as surface 
currents. 
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2.2 Equations for an approximate single coating. 

Next, consider the situation when a metamaterial coating only approximates 
this ideal invisibility coating. We show that the existence of the surfac cur- 
rents for the ideal cloak causes a blow up of the fields as the approximation 
tends to the singular ideal material e, \i. 

To this end, we modify the construction described in the previous sec- 
tion, still dealing with a cloaking structure of the single coating type. More 
precisely, for 1 < R < 2, consider an infinite cylinder in M 3 given, in cylin- 
drical coordinates, by N 2 = {r < R}. On we choose the metric to 
be Euclidian, so that the corresponding permittivity and permeability are 
homogeneous and isotropic. In M. 3 \ N^, we take the metric g and the cor- 
responding permittivity and permeability e and Jl to be the single coating 
metric considered in [6], and the previous section, truncated by being re- 
stricted to N^. Thus, we start from the materials e and Jl corresponding to 
the single coating metric g outside and and replace the metric with the 
Euclidian metric in N^. Then the anisotropy ratio, 

L R := sup ( max =0((R~ l)" 2 ) - oo, 

as the approximate cloaking construction approaches the ideal, that is, R — > 
1+ 

Next, we consider the wave propagation phenomena that arise as the 
approximate cloaking construction approaches the ideal. 

3 Analysis of solutions 

Assume that k is not a Neumann eigenvalue for the Euclidian Laplacian in 
the 2— dimensional disk {r < 1}; as will be seen later, this is equivalent with 
the condition (J )'(k) ^0. For 1 < R < 2 fixed, let 

N = {r > 2}, 

N* = {R<r< 2}, and 

N? = {r< R}, 

so that the Euclidean space R 3 is the union iV = iVo U N x U N 2 . Let Y, R = 
{r = R} be the (appoximate) cloaking surface and v = d r be its Euclidean 
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normal vector on both sides, T, R . To define the approximate cloaking material 
parameters e R and J1 R , introduce, as in [6], an auxiliary space M R , and 7j = lj 
the Riemannian metric corresponding to this construction. M R is obtained 
by taking the disjoint union of three components, 

M = {r > 2}, 
M R = {p<r < 2}, 
K = {r< p}, 

where p = 2(R — 1). 

The domain M UM R is a subdomain of M. 3 , as is the cylinder M R . Define, 
as in [6], an abstract manifold M R by gluing points (p, 9, z) of the boundary 
of M U Mf with points (R, 9, z) of the boundary of Mf. Equip M R with 
the Euclidian metric g and the corresponding homogeneous, isotropic per- 
mittivity and permeability, e, /i. With respect to the cylindrical coordinates, 
(r, 9, z), we have 

/ 1 \ / r \ 

9 = M?, fc =i = r 2 , e = /i = r- 1 . (2) 
\ 1 / \ r / 

Next, introduce a Lipschitz-diffeomorphism F R : M fi — >• R 3 , which in cylin- 
drical coordinates is given by 



pR . 


;M - 


>N , 


r \Mo - 


id, 


pR . 


: M R - 


^N R , 


F \ m r{ 


r,9,z) = (r/2 + l,9,z 


pR 


;M R - 


^N R 




= id. 



We define the metric g R on R 3 by the formula g R = (F R )*g, that is, 

p,q=l 

where \g 3 ] = [fi^fe] -1 and we use that <? P9 = Suppressing for the time be- 
ing the superscript R , the permittivity and permeability, e, p, corresponding 
to the metric g are then given (see, e.g., [H [2]) by 

e = Jl=\det(g jk )\ 1 / 2 g^. 
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Then the metric g and permittivity and permeability, e, // are still given by 
formula l^j on N and N 2 . On Ni they are 

4 
4(r-l) 2 
1 

/ (r - 1) 
e = /I= (r- l)- 1 

\ 4(r-l) 

In the following, we consider TE-polarized electromagnetic waves. This 
means that, written componentwise with respect to either coordinate sys- 
tem as 
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E = (Ei, E 2 , E 3 ) = (E r , Eg, E z ) 

with 

E r = Ei cos(#) + E 2 sin(0), Eg = r (-Ek sin(0) + E 2 cos(0)) , E z = E 

the electric field has a nonzero component only in the z-direction, 

Ei = E 2 = E r = Eg = 0, E 3 (x) = E 3 (r, 6). 
We denote E 3 = u. Then 

H = -r/2 _1 (V x E) = x Vu). 

IK 

We note that tt satisfies the (scalar) Helmholtz equation, 

(Ag + k 2 )u = on M 3 
where is the Laplace-Beltrami operator corresponding to the metric g. 

3.1 Scattering problem 

We consider an incoming TE polarized plane wave. In N such a wave has 
the form 

E in (r,6,z) = e tkrcose = ( J (kr) + f^2i n J n (kr) C os(n6)\ e z , 
1 



H in (r,6,z) = -^Vxfi 

IK 



n=l 

-u 



or, in terms u in = E 3jin , 

oo 

J (kr) + 2i n Jn(kr) cos(n0) 



Ui,_ 

n=l 



Here fi = e = 1. We look for the solution of the scattering problem, 

V x E = ikB, V x H = -ikD, (3) 
D = eE, B = JmH 

on M 3 , where e = e R , Jl = J1 R so that E = E R , etc. Supressing again 
the index R , E = E in + E sc , H = H in + H sc , and E sc and H sc satisfy the 
Silver-Muller radiation condition [2B]. Analysis of cylindrical cloaking using 
Fourier-Bessel series is also in [20] . 

We recall that JR 3 = N U N 1 UN 2 - In the domain N = {r> 2}, one has 

E sc (r.8,z) = (jTc n HW(kr) C os(n6))e z , 

\n=0 / 

H sc (r,9,z) = ^(Vxi), 



Use 

n=0 



J^CnH^ikr) cos(n9) 



Now use the change of coordinates, F : M — > N to define the pulled back 
fields on M, 

zp p* p rr _ p* rr 

"in "in; "in "in; 

"sc "SC> "SC ± ±± SC, 

E = F*E, H = F*H. 

In the coordinates (V, 0', z') = F -1 (r, 0, z),6 f = 6,z! = z on M = F _1 (iVo), 

Vx£ = V x H = -ikD, 

D = e E, B = fj, Q H. 

In Mi, i.e. for r' > p, 
E(r',6',z') = ( Jo(fcr') + c H^(kr') +f^{2i n J n {kr') + c n H^(kr')) cos(n9')) e z 

\ n=0 / 
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H(r',9',z') = ^(Vx£), 

u(r', 9') = E 3 (r', 9') = J (kr') + c H^\kr') + 

oo 

+ ^(22 n J n (£x') + c n H^\kr')) cos(n0') 

n=0 

E r (r, 9, z) = 2E r (2(r - 1), 9,z), ^ E e (r, 9, z) = E e (2(r - 1), 6, z), 

E t (r,e,z)=E g (2(r-l),9,z). 

In N 2 , i.e. for r < R, 

E(r,9,z) = y^a n J n (kr)cos{n9)\e Z} (4) 

\n=0 / 

oo 

u(r,9) = a n J n {kr) cos(nfl), 

n=0 

ff(r,0,z) = ^(VxS). 

As = the fields E and if and the potential u are also given by (J4]) 

in M 2 . 2 

On S R = dN^, using the standard transmission conditions for the electric 
and magnetic fields, that ensure distributional solutions, are 







= E e \ 




E z \t, r + 


= E Z \-£ R -] 


H e 




= H e 




H z \t, r + 


= H z \z R - 


b r 




= D r 




B t \t, r + 





we get the following transmission conditions for u, 

u\z R + = u\x R -, 
4(R - 1) d r u\T, R+ = Rd r u\Y, R -- 



These correspond to conditions on dM 2 = <9(M U M 1 ), 

pd r u\ r=p + = Rd r u\ r= Fi-, 
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that give equations for c n and a n . 

Let us start with n = 0, which is of particular interest. The above 
conditions yield for ao and cq the equations 

a J (kR) = J (kp) + c H^\kp), 

a Rk(J )'(kR) = pk(J o y(kp) + c pk(Hl ) 1) y(kp) 

that yield, when (J )'(k) ^ 0, that 

p(j y(kp)j (kR)-Rj (kp)(j y(kR) . i f m . 

p{H^y{kp)J (kR)-RH^\kp)(J y{kR) MM 

a = kp{J y{kp)H$\kp) - kpJ (kp)(H^)'(kp) = 2vr + 

"' p(H^y(kp)J {kR)-RH^\kp){J y{kR) ( Jo)'(^) MM 

where we use the asymptotics of Bessel functions near 0, see [23, pp. 360- 
361]. Here, o(l) means that the quantity goes to zero as i? — > 1 + , p — > + ,. 
Similarly, a n and c n satisfy the equations, 

a n J n (kR) = J n (kp) + c n H^\kp), 
a n Rk{J n )\kR) = pk(J n y(kp) + c nP k(H^y(kp), 

that yield, for the generic k, that 

c (jR) = P (J n y(kp)UkR) - RJ n (kp)(J n y(kR) = 
P {H^)'{k P )J n {kR)-RH^{k P ){J n )'{kR) 

a (R) = k P (j n y(kp)Hi l) (kp) - kpj n {k P ){H^y{kp) = 

p(H^y(kp)J n (kR)~RHL 1 \kp)(J n y(kR) 

This implies that the scattered fields (far-field patterns) E sc , H sc in N U N R 
and the transmitted fields E, H in N R , which, as we recall, depend on R, go 
to zero as the approximate cloaking construction tends to the ideal material 
parameters, i.e. R — > 1 + . A similar result was obtained in [20] . 

Next, we consider the behavior of the fields E R , H R , D R , B R near = 
dN R . Supressing again the superscript R , we write the electric and magnetic 
fields as 

oo 

E(r,9,z) = ^E n (r,9,z), where E n {r, 9, z) = f n (kr) cos(n9)e z ; (6) 

n=0 

H(r, 9,z) = V H n {r, 9, z), where H n {r, 9, z) = —p' 1 (v x E n {r, 9, z)) , 



n=0 
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with similar notations for the scattered and incoming fields, E™ c , H™ c , etc. 
On M, the decomposition gives rise to a similar decomposition of E 
and H, which we analyze for each value of n. First, we consider the terms 
corresponding to n = 0. On M U M 1 , at y = F{ 1 (x),y = (r',0',z'), x E 
N U Tvf , 

ElM = J (kr')=O(l), 

E° SC ,M = c (R)H£\kr>) = - 1 -^ (l + o(l)). 

Observe that, since r' > p, E® cz (y) is uniformly bounded for R — > 1 + . With 
the magnetic field H° having a non-zero component only in 6, one has 

Hl, 6 {y) = ir' (Jo)' (kr') = 0((r'n 

HIM = ir'c (R)(H^)' (kr') = ^-^(1 + 0(1)). 

On Mf , 

Ofl) 



E° z (y) = a (R)J (kr') 



H° e (y) = ia Q (R)r' (J )' (kr>) 



In (kp) ' 

0(1) 



In (kp) 



Returning to iV and again using the transformation rules for E and H, we 
see that E°, H° are uniformly bounded, with respect to R, in N 1 U N 2 . 

Now consider the magnetic flux density, B = pH which has a similar 
decomposition. In particular, on N^, one has 

BlM 9) = pHlM 6) = pHl e (2(r - 1),6) = 0(r - 1), (7) 
B%(r, 6) = pH° sc /r, 6) = pH^ (2(r - 1), 9) = * (1 + o(l)). 



Pointwise, on N 2 , 



n°(rf))- 0(r) 
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tending to when R — > 1 + . However, to see how Bg behaves as a distribu- 
tion as R — ► 1 + , observe that implies that J^ 2 Bg(r,9) dr is uniformly 
bounded as i? — > 1 + , while for any < k < |, 

1+K * f ^ 



£°(r, 9)dr = - / -dt + o(l) -> — when i2 = p/2 + 1 -> 1+. 

i-« k J p/2 (log p)t k 



This implies that 



lirn fl° = ^ s + B° b>(h 



in the sense of distributions, where 8% is the delta-function of the cylinder 
E = {r = 1} and is a bounded function. 

At last, consider D° which has only the z— component different from 0. 
In Nf, 



Dl >z (r, 9) = eEl z (r, 9) = (r - l)E^ z {2{r - 1), 0) = 0{r - 1); 
5° (r, 9) = eE°(r, 9) = (r - l)E° sc J2(r - 1), 5) 



o , a , , , nt _ ^ ^ 0((r-l)ln(r-l)) 



In (fcp) 



while in iV<f , 



D°(r, 9) = eE° z {r, 9) = rE° z {r, 9) ° l ' 



In (kp) 

Thus, when R — > 1 + , Df n2 has a uniform limit in NoU N±, and .D° cz , -D° 
uniformly tend to in jV U iVi , N 2 , respectively. 

For n > 1, using (jSJ), we obtain the following asymptotics for E 1 , etc in 
various subdomains of iV: 

In iV/ 2 , where r > R, i.e. 2(r — 1) > p, 



p 2n 



r - i; 

2n 
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/ n2« 

D? n , z = 0((r-ir +1 ), d: c =o' ' 



(r - l) n 

B: n =0((r-lT), b: c =o 



. J / 

p 



r - i; 

2n 

n-l\ 



V _ l)n+l 



^^^((r-ir- 1 ), 5^ = 
As for N R , we have 

H^ = 0(p n ), HZ = 0(p n ); 
B? = 0(p n ), B% = 0(p n ) 

These formulae imply that there is a uniform limit of E n , H n , D n , B n when 
R — > 1 + and, moreover, the scattered fields in iVo U N\ and transmitted 
fields in N 2 tend to 0. These formulae also imply that the series 

oo oo oo oo 

5Z-^™ c ' 5Z-^™ C > X^-^" c ' 

n=l n=l n=l n=l 



in iV U JVi; and 

oo oo oo oo 

j2E n , Y,H n , E^ n ' E 5 "' 



n=l n=l n=l n=l 



in N 2 , all converge to zero, as R — > 1 + , in the sense of distributions, i.e., in 
V'(N, dx). 

Summarizing, we see that, in the sense of distributions, 



lim E R = E b , lim H R = H b , 

R-*l+ R^l+ 



1 

r^T+ " ik 



lim D R — D b — -TJ-Jsurfi Jsurf — 0, 



~ ~ 1 ~ ~ 
lim B R = B b + —K mr f, K SUT} = -5%. 

Here E b , H b , D b , and B b coincide with E in , H in , B in , and B in , correspond- 
ingly, in No U Ni and are equal to in N 2 . Thus, in particular, they satisfy 
equations ([3]) separately in N U Ni and N 2 . 
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Note that, sending an TM-polarized wave, we get that J sur f = — <5s, K sur f = 
0. Moreover, for a general incoming electromagnetic wave, the correspond- 
ing solutions E R , H R , D R , B R tend in L^M 3 ) to E Um , H Um , D lim , B Um which 
satisfy 

V x En m = ikBn m + K sur f, V x Hn m = ikDn m + J SU rfi 
with J surf = 5 e <5 E , K sur f = b h fe. 



4 Numerical results 

We next use the analytic expressions found above to compute the fields when 
a plane wave, with vertically polarized E-field, E in (r, 8, z) = e tkrcose g z having 
wavenumber k = 3, is incident to a cylinder {r < R} that is coated with 
an approximative invisibility cloaking layer located in {R < r < 2}. We 
then numerically simulate the cases where R = 1.01 and R = 1.05. In 
the simulations we have used Fourier series representation to order 6, that 
is, the fields are represented using trigonometric polynomials of degree less 
than or equal to six, ^i n |<6 fn( r )e m9 - In the tables below, we give the real 
parts of the ^-component of the total fields and the scattered 5-field on the 
line {(x, 0,0) : x G [0,3]}, first in the absence of a physical layer inside 
the metamaterial and then when an SHS lining is included. We note that 
in the case of the SHS lining, the fields are as was claimed of [U [7J [H [T7] 
without reference to a lining, namely zero inside the cylinder {r < R}. In 
Figs. 1,2, we see clearly the development of a delta-type distribution on 
the interface when we do not have the SHS lining and the approximative 
cloaking approaches the ideal, i.e., R — > 1 + . Also, we see that far away from 
the coated cylinder in both cases the scattered field goes to zero, but much 
more quickly when the SHS lining is used. 
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Below we give the numerically computed Fourier coefficients of the scat- 
tered waves. When L = 1600, that is, R = 1.05, we have: 



n 


c n with SHS lining c n 


with no SHS lining 





-0.0042 - 0.0644i 


-0.6308 - 0.4826i 


1 


-0.0099 + 0.1407i 


-0.0067 - 0.1154? 


2 


-0.0016 + O.OOOOi 


-0.0008 - O.OOOOi 


3 


0.0000 + O.OOOOi 


0.0000 + O.OOOOi 


4 


0.0000 + O.OOOOi 


0.0000 + O.OOOOi 


5 


0.0000 + O.OOOOi 


0.0000 + O.OOOOi 


6 


0.0000 + O.OOOOi 


0.0000 + O.OOOOi 


(EM 2 ) 


1 / 2 0.1551 


0.8026 


When L = 


= 40, 000, that is, R = 1.01, we have: 


n 


c n with SHS lining c n 


with no SHS lining 





-0.0000 - 0.0028i 


-0.2318 - 0.4220i 


1 


-0.0000 + 0.0057? 


-0.0000 - 0.0049i 


2 


0.0000 + O.OOOOi 


0.0000 + O.OOOOi 


3 


0.0000 + O.OOOOi 


0.0000 + O.OOOOi 


4 


0.0000 + O.OOOOi 


0.0000 + O.OOOOi 


5 


0.0000 + O.OOOOi 


0.0000 + O.OOOOi 


6 


0.0000 + O.OOOOi 


0.0000 + O.OOOOi 



(EM 2 ) 1/2 °- 0063 °- 4815 

The results show that for R close to 1, including the SHS lining strongly 
reduces the farfield of the scattered wave; the approximative invisibility 
cloaking functions much better with such a lining than without, even for 
cloaking passive objects. 

5 Discussion 

5.1 Comparison of results with and without SHS. 

One observes that, without the SHS lining, the 5-field grows as the approx- 
imate single coating tends more closely to the ideal invisibility cloak, i.e., as 
the anisotropy ratio L becomes larger. In both Figs. 1 and 2, the peak near 
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0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



Figure 1: The ^-component of the total I?-field on the line {(x, 0,0) : x G 
[0, 3]}. Blue solid line correspond the field with no physical lining at {r = R}. 
Red dashed line correspond the field with Soft-and-Hard lining on {r = R}. 
At the left figure, the maximal anisotropy ratio is L = 1600 and R = 1.05. At 
the right figure, the maximal anisotropy ratio is L = 40, 000 and R = 1.01. 




Figure 2: The ^-component of the scattered 5-field on the line {(x, 0,0) : 
x G [0,3]}. Blue solid line correspond the field with no physical lining at 
{r = R}. Red dashed line correspond the field with Soft-and-Hard lining on 
{r = R}. At the left figure, the maximal anisotropy ratio is L = 1600 and 
R = 1.05. At the right figure, the maximal anisotropy ratio is L = 40,000 
and R = 1.01. 
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r = 1 without the SHS lining shows quite clearly how the delta-distribution 
in the 5-field develops. 

Note that the value of the anisotropy ratio L is quite large in our sim- 
ulation, but the resulting fields are still not extremely large. So, it is not 
surprising that the non-existence results in [6], predicting the blow up of 
the fields, were not observed in the experiment [S] . However, it seems likely 
to become more significant as cloaking technology develops. Also, the SHS 
boundary lining has the additional benefit in our simulations of making the 
scattered wave smaller outside of the metamaterial construction. Indeed, the 
scattered field when using the SHS boundary lining is less than 2% of the 
scattered field without the lining. Thus, implementation of a lining signifi- 
cantly improves the cloaking effect. 

5.2 Significance of the surface currents J sur f and K sur f 

As R —>■ 1 + , for generic incoming waves the magnetic and electric flux den- 
sities converge to fields that contain delta-function type distributional com- 
ponents supported on the surface E. We phrase this by saying that surface 
currents appear. If the metamaterial construction allows this, then we inter- 
pret this literally. This holds, e.g., if the metamaterials used have components 
near E that approximate a SHS surface, such as strips of PEC and PMC ma- 
terials. Alternatively, if no such currents can appear in the material, the D 
and B fields will blow up as the approximation of the coating material goes 
to the limit 

Effective medium theory for composite materials is proven only when the 
limiting fields are relatively smooth [28J. Such rigorous effective medium 
theory has not yet been established for metamaterials, but the limited work 
so far, e.g., [29J, clearly indicate that this same restriction will hold there as 
well. One can then interpret the blow up of fields as a challenge to the valid- 
ity of the material parameters that have been ascribed to the metamaterilas 
currently employed. Indeed, fields having a blow up are very rapidly chang- 
ing functions near the cloaking surface. Thus making a physical cloaking 
construction that would operate well with such fields would require metama- 
terials whose cell size becomes very small close to the cloaking surface. 

The simplest way to avoid these issues would be to include the SHS lining 
when constructing the cloaking device. 
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5.3 Summary 

We have considered two cases when cloaking an infinite cylinder: 

(1) An infinite cylinder of air or vacuum, is coated with metamaterial in {R < 
r < 2} but has no lining on the interior surface of the metamaterial coating. 
In the limit R — > 1 + , solutions to Maxwell's equations have singular current 
terms K sur f and J SU rf that represent either surface currents or the blow up of 
the D and B fields. A standard assumption in homogenization theory is that 
the length scale, d, of the substructures (or cells) from which a composite 
medium is formed, is much less than the free space wavelentgnt A of the EM 
field [28J. In treatments of homogenization for metamaterials, e.g., [29], it has 
been observed that effective material parameters can often be obtained even 
when d is not greatly less than A. Although not explicitly stated, it is required 
that sampled surface intergrals of E, H, B, and D not vary greatly from point 
to point within a metamaterial cell. The blow up of B that we have shown 
occurs when cloaking wihout an SHS lining thus presents a challenge to the 
effective medium intepretation of the metamaterials employed. 

(2) An infinite cylinder of air or vacuum is coated with metamaterial in 
{R < r < 2} and a SHS-lining on the interior of the cloaking surface. The 
lining can be considered as parallel PEC and PMC strips, that allow surface 
currents in the ^-directions. In this case, when R — > 1 + , the total E and H 
fields at the boundary have very small ^-components, that is, in the limit the 
tangent components of E and H are ^-directional. The non-zero tangential 
boundary values of E and H correspond physically to surface currents, that 
are now allowed because of the SHS lining. Since the surface lining and fields 
are now compatible, the fields do not blow up. In addition, the amplitude of 
the farfield pattern is greatly reduced. 
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